Hysteresis and complexity in the zero-temperature mean-field RFIM: the soft-spin 

version 



M. L. Rosinberg^E and T. Munakata 2 'Q 

1 Laboratoire de Physique Theorique de la Matiere Condensee, CNRS-UMR 7600, 
Universite Pierre et Mane Curie, 4 place Jussieu, 75252 Paris Cedex 05, France 
2 Department of Applied Mathematics and Physics, 
Graduate School of Informatics, Kyoto University, Kyoto 606-8501, Japan 

We study the energy landscape of the soft-spin random field model in the mean-field limit and 
OS \ compute analytically the quenched complexity of the metastable states as a function of their mag- 

netization and energy at a given external magnetic field. The shape of the domain within which 
the complexity is positive (and the number of typical metastable states grows exponentially with 
system size) changes with the amount of disorder and becomes non-convex and disconnected at low 
(— ( , disorder. As a consequence, phase transitions occur both at equilibrium and out of equilibrium along 

the saturation hysteresis loop. We focus on the zero complexity curve in the field-magnetization 
plane and its relationship with the hysteresis loop. We also study the response of the system when 
, the magnetization is externally controlled instead of the magnetic field. The main features of the 
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model that should survive in finite dimensions are discussed. 
PACS numbers: 64.60.av, 75.60.Ej, 75.10.Nr 

I. INTRODUCTION 



The response of systems with quenched disorder to an applied field or force is typically hysteretic and noisy. A 
well-known example is the Barkhausen noise in ferromagnetic materials that results from the intermittent motion of 
, magnetic domain walls in response to a change in the external magnetic field [l| . This behavior, which is observed 
in a wide variety of physical systems, from superconductors Q to martensitic shape- memory alloys Q, can be related 
to the complicated structure of the energy landscape, with a huge number of local minima (or metastable states) 
separated by large energy barriers, which makes thermal fluctuations irrelevant on experimental time scales. A jump 
or avalanche then corresponds to the disappearance of a local minimum in the landscape as the external field is 
j— ( • changed. Much effort has been devoted in recent years to describe the avalanche statistics in various models and to 
explain the origin of the power-law distributions that are ubiquitously observed. Depending on the model and the 
details of the dynamics, scale-invariance results from self-organized criticality or requires fine tuning of the system 
• parameters [3]. 

In the theoretical description, one usually considers a zero-temperature dynamical evolution, which implies that the 
system, prepared in some known initial configuration, visits a deterministic sequence of metastable states (similarly, 
at the thermodynamic equilibrium, the system jumps from one global energy minimum to another while the field is 
changed). Therefore, the stochastic character of the output signal only reflects the intrinsinc disorder in the system 
and the whole information about the noncquilibrium response, for instance at criticality, is encoded in the energy 
landscape. This means that it should be possible in principle (but maybe nontrivial in practice) to describe the 
disorder-averaged properties of the response using a statistical description of the metastable states, that is without 
taking into account the entire history of the system. Of course, this supposes that the nonequilibrium path can 
be characterized unambiguously, for instance by some "extremal" property. This is illustrated in the present work 
where we study a simple mean-field model for which a complete analytical description of the distribution of the 
metastable states can be reached. This model was introduced in Rcf.Q as the starting point for a renormalization- 
group description of the noncquilibrium, zero-temperature random- field Ising model (RFIM), and it may be also the 
basis for a similar study of the energy landscape. 

The zero-temperature Gaussian RFIM with a metastable dynamics is a prototype for a large class of disordered 
systems with avalanchelike behavior [f|. The main feature in three and higher dimensions is the existence of a nonequi- 
librium critical point that separates two different regimes of avalanches [7]. Above the critical disorder, all avalanches 
(i.e. Barkhausen jumps) are of microscopic size and the saturation hysteresis loop, obtained by cycling the magnetic 
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field adiabatically from large negative to large positive values and back, is smooth in the thermodynamic limit. Below 
the critical disorder, there is a macroscopic avalanche at a certain field and the hysteresis loop is discontinuous. In 
the following, we consider a soft-spin version of the RFIM where a spin can take any value between — oo and +00, 
and we study the mean-field limit where every spin interacts equally with every other spin. The model can then be 
viewed as a collection of identical bistable (or Preisach) units interacting via a mean-field term. This type of model 
is very popular in the hysteresis community (see e.g. Ref. 12]) and its properties have been studied in great detail. 
However, as far as we know, the issues that we want discuss here have not been addressed. 

The relationship between the hysteresis loop in the zero-temperature RFIM and the distribution of the metastable 
states in the field-magnetization plane was investigated analytically and numerically in a series of recent papers [1, 
B Ell EH ■ The information about the distribution of the metastable states in the H — m plane is encoded in the 
"quenched" complexity Sq(to, H) which is the logarithm (divided by the number of spins) of the typical number of 
metastable states at the field H with a given magnetization per spin m. The number of metastable states grows 
exponentially with the system size when Eq(to, H) is positive, and it appears that the curve £q(jto, H) = exactly 
coincides with the saturation hysteresis loop in the large-disorder regime || 0, The situation at low disorder is 
more complicated but more interesting, and it has been conjectured that the boundary of the domain of existence of 
the metastable states coincides with the hysteresis loop obtained by externally controlling the magnetization instead 
of the magnetic field. This is an issue which is also relevant to some experimental situations (l2l. [T3j] . Our aim in this 
work is to check analytically this whole scenario in the mean-field model. We shall also consider energetics aspects 
which play an important role at low disorder where phase transitions occur, depending on the energy of the metastable 
states. This information is encoded in the magnetization and energy dependent complexity Eg(m, e, H). 

The outline of the paper is as follows. In Sec. II, we present the model and briefly review the hysteretic behavior 
discussed in Ref.[f|. In Sec. Ill, we compute the complexity, first as a function of magnetization only and then as 
a function of magnetization and energy. We focus on the low-disorder regime where phase transitions occur, both 
in equilibrium (i.e. in the ground state) and out of equilibrium along the hysteresis loop. In Sec. IV, we compute 
the nonequilibrium response obtained by controlling the magnetization and discuss its relationship with the usual 
hysteresis loop and the complexity. Section V concludes with a brief discussion. 



II. MODEL AND HYSTERESIS LOOP 



We consider a collection of N soft spins interacting via the Hamitonian 

n = -jxY, s i s i-Y,( H+hi > i + Y, v ( si ) c 1 ) 

where J > is a ferromagnetic coupling, H is an external uniform field, and {hi} is a set of quenched random fields 
drawn independently from a probability distribution V(h) (in practice, this will be a Gaussian distribution with zero 
mean and standard deviation A). V(si) is a double-well potential that mimics the two states of the hard-spin model. 
As in Ref.[5|, we will choose 

V(s) = \-& + ]t f rS< ° (2) 
I f( s - 1) for s > . 

with k > J (this condition ensures that the magnetization remains finite at any field £f[H). The "metastable" states 
are the local minima of the Hamiltonian where each spin satisfies dTi/dsi = 0, i.e. 

H + Jm + h l 

Si - sign(si) = (3) 

where m — (1/N) J2i s i 1S the magnetization per spin. Clearly, this equation has no negative (resp. positive) solution 
for hi > —H — Jm + k (resp. h < —H — Jm — k). On the other hand, s; can be either positive or negative for 
—H — Jm — k < hi < —H — Jm + k. As a consequence, there is a field interval where the number of metastable states 
with magnetization m grows exponentially with N (and the corresponding complexity is positive). This contrasts 
with the mean-field Ising model studied in Ref. 10] where this number is finite (and actually very small) in the 
thermodynamic limit. 

At T = 0, when increasing the external field adiabatically from —00, a spin remains negative as long as Eq. (fTJ 
admits a solution with < 0. When this is no more possible, the spin moves to the "up" potential well (s, > 0) and 
stays at the bottom of the new well. This move is equivalent to a spin flip in the hard-spin model. Therefore, for given 
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field H and given magnetization m, all spins with hi < —H — Jm+k are negative and all spins with hi > —H — Jm + k 
are positive. The self-consistency condition for the average magnetization then reads 



m T (fi) 
which yields Q 



- H— Jmy+k 



dhV{h) [(H + J m] +h)/k- 1] 



-H-Jm r +k 



dhV{h) [(if + Jm T + h)/k+ 1] 



m T (ff) 



H + k 



2k 



k-J k- J 



- H — Jm v -\-k 



dhV(h) . 



Similarly, when decreasing H from +oo, the magnetization is solution of the self-consistent equation 

H + k 2k r -n-j mi -h 



mi (H) 



k-J k-J 



dhV(h) 



(4) 



(5) 



(6) 



Depending on the disorder strength, these equations admit one or three solutions in a certain interval of H and 
the resulting hysteresis loop is then smooth or discontinous in the thermodynamic limit (the magnetization can only 
vary monotonously with the field and the "unstable" solutions of Eqs. ([5]) or ^ that correspond to a negative 
susceptibility dm/dH have no meaning in this context). Unlike the hard-spin mean- field model, hysteresis is always 
present. Specifically, for a Gaussian distribution V(h), the loop is smooth for A > A c and discontinuous for A < A c , 
with A c /J = \j2j-K kj(k — J)0j. This is illustrated in Fig. 1. (Hereafter, without loss of generality, we shall always 
choose k — 3 in the numerical calculations and take J = 1 as the energy unit.) 





FIG. 1: Mean-field hysteresis loop for the soft-spin version of the Gaussian RFIM with k = 3 and (a) A = 3, (b) A = 0.8 
(A c ps 1.197). In (b), the magnetization has a jump at H — +H c (0-8) ~ ±2.152. The quenched complexity T,Q(m,H) is 
positive inside the shaded area and vanishes on the boundary. The red curve represents the ground- state magnetization. (J is 
taken as the energy unit in all figures). 

In the low-disorder regime, the jump in the ascending branch occurs at a field H = H C (A) where the slope dm^ / dH 
diverges (H C (A) is solution of the implicit equation 2kJP(H c + Jm^(H c ) — k) = k — J). Note that A c is also 
the critical disorder at which a discontinuity first occurs in the equilibrium magnetization curve. The ground-state 
magnetization rriQs{H) is readily obtained by noting that the lowest energy at the field H is obtained when all spins 
with hi < —H — Jm are negative and all spins with hi > —H — Jm are positive. This yields the self-consistent 
equation 



rriGs(H) 



H + k 

k-J 



2k 



J 



-H-Jmcs 



dhV{h) 



(7) 



which also admits three solutions in a certain field interval for A < A c . By definition, the ground state is the state 
with the lowest energy, whence the discontinuity in mas(H) at H — (and the intermediate branch is unstable). 
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III. COMPLEXITY 



Whereas the computation of the quenched complexity is a difficult task when a finite number of spins interact (see 
e.g. [1, El), ^ becomes trivial in the mean-field limit. The results, however, nicely illustrate the qualitative behavior 
of the complexity and its relationship to the hysteresis loop (there are also some special features of the mean-field 
model which will be pointed out in the following). 



A. Complexity as a function of field and magnetization 

We first consider the typical number of metastable states with given magnetization per spin m at the field H, 
irrespective of their energy. The associated complexity £g(m, H) is the quantity that is immediately related to the 
hysteresis loop, as discussed in Refs.§, H, [ll|. By definition, 



E Q (m,H) = km -ln7V(m,F) = linj i - lim - N(m, H) n - 1 



(8) 



where W(ra, H ) is the number of metastable states with magnetization m at field H and the overbar denotes the 
average over disorder realizations (as usual, the order of the limits N — > oo and n — > has been inverted). In the 
present case Af(m, H) is given by 



M{m, H) = J\{ dsi J[ 5( Si - *< - Nm) 



(9) 



where s* is a solution of Eq. (J3j) (from a technical point of view, the simplifying feature of the double-well potential 
V(s) described by Eq. (J^J) is that there is no solution of dH/dsi — 0, i.e. no stationary point of the Hamiltonian that 
is a local maximum) . Introducing the integral representation of the second delta function in Eq. ([9]) and averaging 
over disorder gives 



(10) 



where the replica label a runs from 1 to n and g a is a Lagrange multiplier associated to the constraint sf = Nm 
in replica a. Using Eq. ©, we then obtain 



Af{m,H) n 



J]_dg a 



,AT[A({ 9a })-m 



(11) 



where 



A({ 5a })=ln{ / dhV(h)eZ° B °<> H+Jm+h)/k Y[[e 9 °Q(H + Jm + h + k) + e - 9 °0{-H- Jm- h + k)]} (12) 



and 0(x) is the Hcavisidc function. In the large- N limit the integrals in Eq. (fTT|) are dominated by the stationnary 
points {g*} that maximize A({g a }) — m^2 a g a and are solutions of 



m 



dH{g* a }) , 
dg* lm 

e -A({ 9 :}) f dhV{h)e^3'AH+J m +h)/k 



H + Jm + h + k „»~,„ T , 
e 9 -Q{H + Jm + h + k) 



H + Jm + h-k _ a * 
k 6 



@(-H - Jm-h + k) 



Y[\e 9 '*e(H + Jm + h + k) + e~ 9 ^e(-H - Jm-h + k) . (13) 



Since all replicas are equivalent, we can set g* = g* and take the limit n — > straight away (from now on the 
superscript * will be dropped in order to simplify the notations). This yields 



f „ n ( H + Jm + h + k)e 9 Q(H + Jm + h + k) + (H + Jm + h - k) e - 9 Q(-H -Jm-h + k) 
*™ = y*MW e9 e(H + Jm + h + k) + e- 9 e(-H- Jm-h + k) (W) 
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whence 



H + k 2k f-H-Jm-k 



k-J k-J< < lhnh) 
This finally allows us to express g as a function of m and H 



7 p — H— Jm+fc 

- -[tanh(ff)-l] / dW(fc) 



tanh 5(771, H ) = 1 



(k-J)m-(H + k)+2kf_* Jm k dhV(h) 



(15) 



(16) 



an equation which may have no solution for certain values of H and m since |tanh(g)| < 1. The corresponding 
quenched complexity is obtained from Eqs. © and (fTl]) as 



£ Q (m, H) = A {1) (g(m, H), H) - mg(m, H) 

where 

AW(g,g)=lim A(g ' g) 

ri— >0 77 

= lim - In i 1 + 77 / d/i-pf/i) 
rwo n [ J 



(17) 



g H+J ™ + h + In [e 9 9(iJ + Jm + h + k) + e - 9 S(-H -Jm-h + k)] 



fi- 



ll + Jm 
~k 



+ J dhV(h) In [e g O(H + Jm + h + k) + e - g G(-H - Jm -h+ k)] 



9- 



H + Jm + k 



23 



—H—Jm—k p — H— Jrn-\-k 

dhT(h) + [ln(2cosh(g)) - g] dhV(h) 

00 J—H—Jm — k 



Using Eq. (fT5"]) . we finally obtain the simple result 



/-H-Jm+k 
dhV{h) 
-H-Jm-k 



(18) 



(19) 



where 5(777, H) is solution of Eq. fl 

3 




g(m,H) 

FIG. 2: g(m, H) as a function of m for A = 0.8 and different values of H. For clarity, the magnetization is put on the vertical 
axis. 

The main feature is that Y>Q(m,H) — > as g — > ^00, and one can see immediately from Eq. (|16|) that Eqs. 
([5]) and © that describe the ascending and descending branches of the saturation hysteresis loop, respectively, are 
recovered in these limits. Therefore, the complexity is not only positive in some region inside the hysteresis loop 
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(as a consequence of the no-passing rule[14[), but it is positive everywhere inside the loop for A > A c and exactly 
vanishes along the two branches m^(H) and m^(H). For A < A c , this is only true before and after the jump in the 
magnetization. In between, the curve £q(771, H) = is reentrant, as shown in Fig. 1(b). This part of the curve is 
also described by the "unstable" solutions of Eqs.© and ©, but this is a peculiarity of the mean- field limit [la]. 

The behavior of g{m, H) and Sg(m, H) for A < A c as a function of m and different values of H (on the positive 
side) are shown in Figs. 2 and 3, respectively (in Fig. 2 we actually plot m as a function of g since this representation 
is more intuitive). One can see in Fig. 3 how the shape of T,Q(m, H) is changing as H increases [l8j. In particular, 
for H/J — 2, there are no metastable states with magnetization 1 < m < 2.35. The fact the magnetization of the 
metastable states cannot assume any value is at the origin of the jump in the hysteresis loop displayed in Fig. 1. 




FIG. 3: YiQ(m,H) as a function of m for A = 0.8 and different values of H. Note that the magnetization that corresponds to 
the maximum of EQ(m, H) (i.e. the typical magnetization of the metastable states at the field H) is not a monotonous function 
of H. For H — 2.149, the two maxima of Eq(771, H) have the same height and the typical magnetization jumps discontinuously. 



The typical (i.e. most probable) magnetization at the field H is the magnetization for which Eq(?ti, H) reaches its 
maximum. It turns out that this cannot be simply obtained by setting g = in Eqs. (TIB"]) and (|18| because of the 
explicit dependence of A on m (the only exception is for H = 0, by symmetry). Indeed, whereas dA^(g, H)/dg = m 
(see Eq. CE3D), 

dZ Q (m,H) _ dA^(g,H) dg | dA^\g,H) _ ^ 8g_ 
dm dg dm dm dm 



dm 



(20) 



In other words, £<g(m, H) and AW(j, H) are not not mutually connected by a Legendre transform, which differs from 
the situation in finite-connectivity models [1, [ll|. Eq. (f2"01 then gives 



dE Q (m,ff) 
dm 

In particular, 



k — J 

-g— — + 2gJT{H + Jm + k) + J [In (2 cosh( 5 )) - g] [V{H + Jm + k) - V(H + Jm - ft)] . (21) 



dH Q {m,H) 
dm 



-9 



k-J 



2JV{H + Jm±k) 



for g — > ±oo 



(22) 



so that the curve 9Sg(m, H) /dm = joins the hysteresis loop exactly at the turning points H = ±H C (A) where 
the slopes of m^(H) and mi(H) diverge and the magnetization jumps. However, the actual typical magnetization is 
discontinuous and jumps at a field that is very close but smaller than H c in absolute value. (For the case considered 
in Fig. 1(b), the jump occurs at H w ±2.149, whereas H c (0.8) = ±2.152. This corresponds to the situation where 
the two maxima of Eg(m, H) have the same height.) 
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B. Complexity as a function of field, magnetization, and energy 



We now complete our description of the distribution of the metastable states by computing the complexity at the 
field H for a fixed magnetization m and a fixed energy per spin e. This is achieved by introducing an additional 
Lagrange multiplier /3 which plays the role of the inverse temperature for the metastable configurations. Since the 
calculation proceeds along the same lines as before, we only quote the result, 

p-H-Jm+k 

S Q (m, e, H) = / dhV{h) {In (2cosh[ 5 + (3(H + Jm + h)}) - [g + /3(H + Jm + h)} tanh [5 + (3{H + Jm + h))} 



H—Jm — k 



(23) 



where g and /3 are solutions of the implicit coupled equations 



H + k 



2k 



k-J k-J 



-H-Jm-k 



dhV(h) 



k-J 



-H-.Jm+k 



dhP(h) [tanh [g + (3 (H + Jm + h)] - 1] 



(24) 



-H-Jm-k 



and 



Jm 2 (H + Jm) 2 + A 2 ,__ _ . n 
e = — (H + Jm) + 2 



2 

-H-Jr, 
-H-Jm-k 



2k 



-H—Jm—k 



dhV(h)(H + Jm + h) 



dhP(h)(H + Jm + h) [tanh [g + (3 (H + Jm + h)] - 1] 



(25) 



Conversely, we may consider that the complexity, the magnetization, and the energy are parametrized by g and (3. 
One recovers the equations of the preceding section when /3 — 0. This corresponds to the maximum of Sg(m, e, H) 
at fixed m and H. A randomly chosen metastable state with magnetization m{g, [3 = 0) will thus have an energy 
e{g,(3 = 0) with probability 1 in the thermodynamic limit. 

From these equation we can now study the domain of existence of the metastable states as a function of H , m, and 
e, and determine the minimum (resp. maximum) value of the energy for fixed H and m below (resp. beyond) which 
there are no typical metastable states and Eg(m,e, H) — (or, alternatively, the minimum and maximum values of 
m at a given energy level e). Eq. (|23|) tells us that T,Q(m, e, H) goes to zero when g and (3 go to ±oo, but the ratio 
r — gj ' (3 must be kept finite so that m and e can be varied continuously by controlling r. One finds from Eq. (I24|) or 
()25p that |7'| < k, and depending on whether (3 — > +oo or — oo these equations become 



H 



2k 



-H — Jm—r 



dhV{h) 



(26) 



Jm 2 (H + Jm) 2 



2k 



(H + Jm)+2 J 



-H—Jm—r 



dhV(h)(H + Jm + h) 



(27) 



H + k 



2k 



k-J k- J 1 



-H-Jm-k 



dhP(h) 



-H-Jm+k 



— H—Jm—r 



dhV(h)] 



(28) 



Jm 2 (H + Jm) 



2k 



-H-Jm-k 



(H + Jm) 



dhV{h)(H + Jm + h) 



-H-Jm+k 



dhV(h)(H + Jm + h) 



— H—Jm — r 



(29) 



Both the hysteresis loop and the ground state are recovered for special values of r. The ascending (resp. descending) 
branch of the hysteresis loop, m^(H) (resp. mi(H)), is recovered from Eq. l|2"6"l) for r = —k (resp. +k) or from Eq. 
([28)1 for r = +k (resp. —k). The equation for the ground-state magnetization, Eq. (JT)), is recovered from Eq. (|26[) for 
r = 0, whereas Eq. (I28| for r = gives the magnetization of the metastable states with the maximum energy. 

In the large-disorder regime A > A c , Eqs. (|26p and (|28J) . considered as implicit equations for m, have a unique 
solution for all values of H and r and the generic behavior of the curve Eg(m, e, H) = in the plane m — e is quite 
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FIG. 4: The curve Eg(m,e, H) — in the magnetization-energy plane for A = 3 and H = 0. The parts of the curve that 
correspond to /3 — > +oo and /3 — > — oo are drawn in black and red, respectively, and the arrows indicate the sense of variation 
of the parameter r. 




m m m 

FIG. 5: Same as Fig. 4 for A = 0.8 and (a) H = 0, (b) H = 1, (c) H = 2. For H = 0, the typical magnetization of the 
metastable states at a given energy, i.e. the magnetization at which the complexity T.Q(m, e, H = 0) is maximum, is plotted 
as a function of energy (blue dashed line). The typical magnetization is non-zero for e < e* ~ —0.533. For H = 2, the domain 
within which Eg(m,e, if) is positive is disconnected (the inset is a blow-up of the little domain on the right). 



simple. This is illustrated in Fig. 4 for A = 3 and H = (we also show in the figure how the parameter r varies along 
the curve). More generally, the behavior of Xq(to, e, H) for A > A c is rather featureless and not worthy of special 
comment. Some typical results in zero external field are shown in Figs. 7(a) and 8(a) below. 

Much more interesting is the behavior in the low-disorder regime A < A c . In this case, Eqs. Prj|) and (J25J) may 
have three solutions in a certain range of H and r and the domain within which Eg(m, e,H) is positive is no more 
convex and may even break into disconnected regions, as shown in Fig. 5. This gives rise to phase transitions. 

Let us first consider the situation in zero external field. A remarkable feature is that the typical metastable 
configurations have a non-zero magnetization at low energies (and hence low temperature). This is shown in Fig. 6 
where we plot the complexity as a function of m for A = 0.8 and three different values of the energy. While the 
complexity has a single maximum at to = for e = 0, it exhibits two maxima at symmetric values of the magnetization 
at lower energies. In this case, the spin configurations that dominate the distribution of the metastable states in the 
thermodynamic limit are magnetized. This occurs below a certain energy e*(A) (or, equivalently, for f3 larger than 
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FIG. 6: Complexity Eg(m, e, H) as a function of magnetization for A = 0.8, H = and three energy values, e = 0, e = —0.7, 
and e = —0.8. At low energies, the complexity has a double-peak structure and the typical metastable configurations are 
magnetized. Eventually, there is a range of magnetizations to which no metastable configurations can be associated. 




FIG. 7: Complexity of the metastable states as a function of their energy e (and irrespective of their magnetization) for H — 0: 
(a) A = 3, (b) A = 0.8. In (b) the typical magnetization of the metastable states is zero for e > e* and non-zero (red curve) 
for e < e* , with e* ~ —0.533. For e < e* , the complexity of the metastable states with m = (dashed curve) is smaller. 



some threshold j3* (A)), as shown in Fig. 5(a). We notice that a similar result was found numerically for the low-energy 
metastable configurations in the 3d RFIMpJl, which shows that this is not an artifact of the mean- field model (the 
same feature also occurs for the pure Ising ferromagnet on random regular graphs (20l l2l| ) . In contrast, the typical 
magnetization of the metastable states is zero at all energies for A > A c . 

When H = 0, the extrema of the complexity and the corresponding magnetizations and energies can be readily 
obtained by setting g — in Eqs. (I23[l - (f55]) because of the symmetry H <-» —H, m <-» — m (this is not true in general, 
as discussed previously). The maximal complexity Sg(e, H = 0) = max m Sq(tti, e, H = 0), i.e. the complexity of the 
metastable states irrespective of their magnetization, is plotted in Fig. 7 as a function of energy, and the corresponding 
typical energy is plotted in Fig. 8 as a function of the inverse temperature f3 (in Figs. 7(b) and 8(b), we also show 
the branches that correspond to the states with zero magnetization which have a smaller complexity). 

The transition from configurations with zero magnetization at high energies to magnetized configurations at low 
energies is a second-order phase transition. The inverse temperature /3*(A) is thus obtained by expanding Eq. (|24|) 
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FIG. 8: Typical energy of the metastable states against inverse temperature for H = 0: (a) A = 3, (b) A = 0.8. In (b), the 
kink in the energy at /3 = j3* « 0.953 corresponds to a second-order phase transition: for /3 > (3* the typical metastable states 
are magnetized (the dashed curve indicates the energy of the metastable states with m = 0). 



(with H = and g = 0) around m = and cancelling the term proportional to m. This yields the equation 

- V(k) + [ dhV'(h) tanh(/3* /i) = (30) 
2kJ J 

where V'(h) = dV(h)/dh. The corresponding energy e*(A) is obtained by setting H = 0, m = 0, g = 0, and 
j3 = /3*(A) in Eq. (|25p. Eq. (j3"0|) thus defines a second-order line in the disorder-energy and disorder-temperature 
planes separating a paramagnetic and a ferromagnetic phase, as shown in Figs. 9 (a) and 9(b). This extends the 
equilibrium phase diagram of the model to metastable states of low energy. 

1| i | , | i | i I | , | , | i | , 1 2 




FIG. 9: Phase diagram of the transition to ferromagnetic metastable states in zero external field. In (a), the energy e* at the 
transition is plotted against the disorder A (black solid line); we also show the minimum (red dashed line) and maximum (blue 
dotted-dashed line) allowed values of the energy. In (b), the transition "temperature" 1/13* is plotted against the disorder. 
At low energy (or low temperature) and low disorder, the typical metastable states are magnetized. Otherwise, their typical 
magnetization is zero. 

In non-zero field, the most salient feature is that the domain within which £g(m, e, H) is positive may be discon- 
nected, as shown in Fig. 5(c). This is of course in agreement with the results of the preceding section, and, as already 
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emphasized, is at the origin of the jump in the hysteresis loop. 



IV. H-DRIVEN AND M-DRIVEN PROTOCOLS 



Another issue that is worth discussing in the present context concerns the influence of the driving mode on the 
nonequilibrium response of th e sy stem. This issue has been discussed in detail in the case of sandpiles or driven 
interfaces in disordered media[22| (for instance when an elastic chain is driven adiabatically at constant force or 
constant velocity on a disorder substrate[23]). Experimentally, when the variable conjugated to the force (or field 
or stress) is externally controlled instead of the force itself, one observes a reentrant hysteresis loop as well as large 
fluctuations of the induced force [12|, [l3|. The magnetization-driven RFIM exhibits this behavior in the low-disorder 
regime [24], and it has been suggested that the resulting hysteresis curve follows the boundary of the domain of 
existence of the typical metastable states in the thermodynamic limit. We can now test this assumption in the mean- 
field model since we have an exact analytical description of the curve £q(to, if) = 0. Note that the problem of finding 
the magnetic field that produces a desired magnetization function is known in the hysteresis literature as the "inverse 
Preisach problem" [26| . 

In the present case, the magnetization-driven (or Af-driven) protocol is defined as follows. As the magne- 
tization to = M/N is slowly varied, the system attempts to minimize (at least locally) its internal energy 
M({ s i}) — s i s j ~ J2i hi s i + J2i V( s i) (there is no externally imposed magnetic field). Each spin Si 

thus satisfies dU/dsi — with the global constraint YliLi s i — Nm. This simple problem of constraint optimization 
is readily solved by introducing the quantity 



C = U -XV^Si- Nm\ (31) 



. i=l 



where A is a Lagrange multiplier that has the meaning of a magnetic field coupled to the extensive variable Yl^—i Sj. 
Minimizing U({si]) with the constraint on the magnetization amounts to solve simultaneously the N + 1 coupled 
equations dC/dsi = and dC/dX = 0, 

, „ l , A + Jm + hi , 
Si {X) - sign( Sl (A)) = - (32) 

with 

N 



^s,:(A)=iVTO. (33) 



Comparing to Eq. ([3]), it is clear that the spin configurations {s*}, solutions of Eq. (f3"2"]l . are mimima of the 
Hamiltonian 7i =U — Hj^i s i f° r the special value of the field if = A*(m) that satisfies the constraint equation, Eq. 
(|33D . The Lagrange multiplier A can be eliminated from Eq. (|32p by summing over i, which gives 



k N 



X = (k — J)m — — ^sign(si) , (34) 

neglecting the contribution (1/iV) ^ hi that vanishes in the thermodynamic limit. Replacing in Eq. (|3"2")) yields 

h 1 N 

St - sign(si) = to + -j- - — sign(sj) . (35) 

(Note that this equation can be simply rewritten as <7i — sign(/j) where /i = Cj + to + hi/k — (1/N)^2aj and 
Oi is the auxiliary Ising variable defined by Ui = sign(si).) The crucial feature in Eq. (|35[) is the presence of the 
antiferromagnetic contribution —(X/N) sign(sj ) which plays the role of an infinite-range demagnetizing field (see 

Ref. [27( for a more detailed discussion of this feature in a slightly different framework) . As is well known, such a term 
may have a dramatic influence on the behavior of magnetic systems II 2jl and models of interface growth in disordered 
media show that it may result in self-organized criticalitvpl [23. l28l. |29|| . 
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There may be of course many solutions to Eq. (|35[) (which simply characterizes all metastable states with magneti- 
zation m) and one needs to define a relaxation dynamics to go from a metastable configuration at m to the "nearest" 
one at m + 5m. The most natural algorithm consists in searching for the fixed point of the iterative map 

s\ n +1) - sign( s |" +1) ) = m + 6m + ht/k - (1/N) ^ sign(sf 5 ) (36) 

j 

applied to all s, in parallel, where sj = s*(m) is the converged value at magnetization m. Since this equation has no 
negative solution when the right-hand-side is larger than 1 and no positive solution when it is smaller than —1, the 
evolution of the spins (and, as a consequence, of the induced field) is discontinuous. Remarkably, it turns out that 
convergence to the fixed point is reached after only one or two iterations. The first iteration gives 

s^ 1 ' — + 5m + [sign(s- ^ 4- 5m) — sign(s| ^)] (37) 

which implies, e.g. when increasing m from a large negative value, that the spins in the range -5m < sf } < become 
positive. However, these spins never become negative again (i.e. there are no "back-flips"). The second iteration then 
yields 

4 2) = a? - (1/iV) £[ S ign( S «) ~ si S n (4 0) )] ( 38 ) 

j 

and there is no need for a third iteration, i.e. s*(m + 5m) = s| 2 . This implies from Eq. (|34|) that the induced field 
at m + 5m along the ascending branch is simply given by 

ff T (m + 5m) = Hf(m) + (k - J)5m - (k/N) ^[sign(s*(™) + Sm) - sign(s*(m))] . (39) 




FIG. 10: Comparison of the nonequilibrium responses obtained by increasing adiabatically the field H (black line) or the 
magnetization m (red line) from — oo to +oo in a single disorder sample of size iV = 50000 for A = 0.8. For clarity, the 
magnetization is on the vertical axis in both cases. The step sizes were 5H — 0.001 and 5m = 0.002, respectively. The ff-driven 
response has a jump whereas the M-driven response is reentrant (the figure (b) is a blow-up in the vicinity of the knee) . The 
blue line represents the curve Eq (m, H) = 0, solution of Eq. ((5} (the discrepancy with the H-driven trajectory before the jump 
is due to finite-size effects). 

A typical example of the nonequilibrium response obtained with this protocol in a single disorder sample in the 
low-disorder regime is shown in Fig. 10 and compared to the result of the if-driven protocol 30] . One can see that 
the M-driven response is reentrant whereas the i?-driven response displays a jump at H « 2.15, in accordance with 
the infinite- N behavior described in section 2 (there are also many small avalanches that can be better seen in Fig. 
10(b)). Moreover, the induced field H(m) follows very closely the curve T,Q(m,H) = that defines the boundary of 
the domain of existence of the metastable states in the thermodynamic limit. They are fluctuations in H(m) but the 
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M-driven trajectory is always inside the //-driven hysteresis loop, as can be seen in Fig. 4(b). Indeed, for in given 



disorder realization, there are no metastable states outside the hysteresis loopjl4j. 

The fluctuations in H(m) decrease with N and we can actually derive from Eq. (|39|) the exact equation of the 
M-driven trajectory in the limit N — » oo, Sm — > 0. Indeed, (1/2 AT) 5ZJsign(s* (to) + Sm) — sign(s*(m))] is just the 
fraction of spins that are negative at magnetization m and that become positive at m + Sm. In the thermodynamic 
limit, the corresponding probability for the random field h is 

-H r (m)-Jm+k dHt(m) 

dhV(h) - [J + n J ]V(Hf(m) + Jm - k) Sm for Sm -»• , (40) 

and inserting this result in Eq. (|39|) we obtain the differential equation 

dH^m) _ k- J- 2kJV(H^{m) + Jm - k) 



dm 1 + 2kV(H^ (m) + Jm - k) 



(41) 



It is easy to see that the solution of this equation is given by Eq. ([5]), considered as an implicit equation for the 
magnetic field as a function of the magnetization, i.e. H^(m^{H)) — H. As shown previously, this equation defines 
the curve Eq(to, H) = in the H — m plane. Therefore, the M-driven hysteresis loop exactly follows the boundary 
of the domain of existence of the typical metastable states in the thermodynamic limit, including along the reentrant 
part in the low-disorder regime, in agreement with the conjecture of Ref. [24j] . Moreover, the //-driven and M-driven 
loops coincide in the large-disorder regime. 



V. SUMMARY AND CONCLUSIONS 



In this work we have studied the soft-spin version of the random-field Ising model in the mean-field limit and com- 
puted analytically the quenched complexity Eg (to, e, H) of the metastable states as a function of their magnetization 
and energy in the presence of an external magnetic field. We have especially focused on the domain T>(m, e, H) within 
which the complexity is positive and the number of metastable states grows exponentially with the system size (and 
the domain T>(m, H) associated to Eq(to, H) = max e Eg (to, e, H) > 0) . The results can be summarized as follows: 

1) At large disorder, the domain T>{m, e, H) is convex. The quenched complexity Eq(to, H) is then strictly positive 
everywhere inside the saturation hysteresis loop and it vanishes along the loop. Moreover, the system responds 
identically when driven by the magnetic field (//-driven protocol) or the magnetization (M-driven protocol). 

2) At low disorder, the domain T>(m, e, H) is non-convex and disconnected. The curve Eg(m, H) = has a reentrant 
part and coincides with the M-driven hysteresis loop. 

3) In zero external field, there is a second-order line in the disorder-energy plane that separates a paramagnetic 
and a ferromagnetic phase: the ground state and the low-lying metastable states are magnetized at low disorder. 

This study provides a pedagogical illustration of the fact that the zero-temperature, nonequilibrium response of 
the system to a slow external driving can be defined unambiguously from a statistical description of the metastable 
states in some appropriate limit. This suggests that the ground state and the metastable states along the hysteresis 
loop can be studied within a unique and purely static theoretical framework, which may offer a possible new starting 
point to investigate the relationship between the critical behavior of the RFIM in and out of equilibrium [3l|. 

To which extent the mean- field picture remains valid when the interaction is short-ranged ? It is of course impossible 
in this case to perform a detailed study of the distribution of the metastable states in energy or magnetization. 
However, previous numerical[l(| and analytical d, |Tl| studies suggest that the essential qualitative features observed 
in the mean-field model survive in finite dimension, in particular the presence of a gap in the magnetization of 
the metastable states at low disorder, gap which is at the origin of the jump along the field-driven hysteresis loop 
(note that in finite connectivity models, this gap is not related to the existence of an "unstable" solution to the 
self-consistent equations describing the hysteresis loop[l6j]: this solution is an artifact of the mean-field limit). There 
is also numerical evidence that the low-energy metastable states of the 3d RFIM in zero field are magnetized at 
low disorder [l9|, suggesting that the phase diagram of the metastable states as a function of energy and disorder is 
qualitatively similar to the one displayed in Fig. 9(a). 

On the other hand, there are also some noticeable new features in finite-connectivity models. Firstly, the critical 
values of the disorder associated to the equilibrium (ground-state) and nonequilibrium (hysteresis) phase transitions 
become distinct |32j]. with A^ vst < AJ? S . Secondly, in the low-disorder regime, the complexity may not be zero along 
the whole boundary of the domain T>(m, H) (this is still putative but supported by analytical calculations on the 
Bethe lattice at the order l/z)[ll|]. Thirdly, the geometry of the lattice may change qualitatively the nonequilibrium 
response of the system obtained by controlling the magnetization. Preliminary simulations on a cubic lattice show a 
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feature that was noticed in previous work[24j and that has been recently discussed in Ref. [27|: when the disorder is 
low enough, the induced field H(m) exhibit strong fluctuations around a certain value that does not depend on 
the magnetization. As a result, the domain V(m, H) has a very peculiar shape in this regime. The behavior around 
Hd actually corresponds to interface configurations in real space, and the fluctuations are similar to those observed 
at a critical dcpinning transition. This is an interesting issue which is worth investigating further. 
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